Is small-world network disordered? 



in 
o 
o 

(N 
> 

o 
o 



I 

C/3 



a 

o 

(N 
> 
(N 

o 

o\ ■ 

o ■ 
^t- : 
o . 
■ 



a : 
o 

o ■ 



X 
S3 



Soumen RojQ and Somendra M. Bhattacharje^j] 

Institute of Physics, Bhubaneswar 751005, India 

Recent renormalization group results predict non self averaging behaviour at criticality for relevant 
disorder. However, we find strong self averaging(SA) behaviour in the critical region of a quenched 
Ising model on an ensemble of small-world networks, despite the relevance of the random bonds at 
the pure critical point. 



I. SELF-AVERAGING 

Very often in physics one encounters situations where 
quenched disorder plays a prominent role. Any physi- 
cal property X of such a disordered system, therefore, 
requires an averaging over all realisations. It would suf- 
fice to have a description in terms of the average [A] av 
where [..] av denotes averaging over realisations ("sam- 
ple averaging") provided the relative variance Rx = 
V x /[X]l v -» for large N, where V x = [X\ v - [X]l v . 
In such a case a single large system is enough to represent 
the ensemble. Such a quantity is called self-averaging. 
Off criticality, if one builds up a large lattice from smaller 
blocks, then thanks to the additivity property of an ex- 
tensive quantity, central limit theorem guarantees that 
Rx ~ A -1 ensuring self-averaging. In contrast, at a 
critical point, because of long range correlations the an- 
swer to the question as to whether X is self-averaging or 
not becomes non-trivial. 

Randomness at a pure critical point is classified as 
relevant or irrelevant if, following the standard defini- 
tion of relevance, it changes the critical behaviour (i.e. 
the critical exponents) of the pure system. Recent renor- 
malization group and numerical studies have shown that 
if randomness or disorder is relevant, then self-averaging 
property is ZosiQ, In particular, R x at the critical 
point approaches a constant as A" — > oo. Such systems 
are called non self-averaging. A serious consequence of 
this is that unlike the self-averaging case, even if the 
critical point is known exactly, statistics in numerical 
simulations cannot be improved by going over to larger 
lattices (large N). 

Let us recollect the definitions of various types of self- 
averaging with the help of the asymptotic size depen- 
dence of a quantity like Rx- If Rx approaches a con- 
stant as N — > oo, the system is non-self-averaging while 
if Rx decays to zero with size, it is self-averaging. Self- 
averaging systems are further classified as strong and 
weak. If the decay is Rx ~ A^ 1 as suggested by the 
central limit theorem, mentioned earlier, the system is 
said to be strongly self-averaging. There is yet another 
class of systems which shows a slower power law decay 
Rx ~ N~ z with < z < 1. Such cases are known as 
weakly self-averaging. The exponent z is determined by 
the known critical exponents of the system. 



The prediction of non-self-averaging nature of critical 
quantities is an extremely significant result coming from 
general renormalization group arguments. This basic re- 
sult of Ref. and the hypothesis of Ref. can be 
summarized as follows. According to finite size scaling, 
when the critical region sets in, the size of the system 
is comparable to the correlation length £ that grows as 
the critical point is approached. The appropriate scal- 
ing variable is N/N c where N c — £ d is the correlation 
volume in d dimensions. At the critical point of a ran- 
dom system, there is an additional source of fluctuation 
from the variation in the transition temperature itself. 
Therefore, instead of the conventional finite size scaling, 
a sample dependent scaled variable is required. A re- 
duced temperature is defined as U = \T — T c (i,N)\/T c 
where T c (i, N) is a pseudo-critical temperature of sample 
i of N sites with T c as the ensemble average of critical 
temperature in the N oo limit . In terms of this tem- 
perature, a critical quantity X is expected to show a 
sample dependent finite size scaling form 



X i (T,N)=NPQ(t i N 1 ^) 



(1) 



where p characterizes the behaviour of [A] av at T C .M 
Thus p = 7/P where v = dv when X is the magnetic 
susceptibility x- The RG approach seems to validate this 
hypothesis especially the absence of any extra anomalous 
dimension in powers of N for Rx ■ Incidentally, this hy- 
pothesis, Eq. ^ excludes rare events of large pure type 
lattices for which pure v should be used. We are not 
considering such cases dominated by these rare events 
(Griffiths' singularity). With this scaling form, the rel- 
ative variance Rx at the critical point or in the critical 
region is given by 



Rx ~ [{5T c )\ v N 2 /\ 



(2) 



where [(<5T c ) 2 ] av is the variance of the pseudo-critical 
temperature. A random system can have several tem- 
perature scales, namely (T C (N) — T c ) and (T — T c ), in 
addition to the shift in the transition temperature itself. 
It is plausible that for a system with relevant disorder 
all these scales behave in the same way so that typical 
fluctuations in the pseudo-critical temperature is set by 



[1] Conventional notations of critical exponents are used: C ~ t~ a , 

X ~ t~ 7 , £ ~ t~ " where C, x an d 5 denote the specific heat, 

'Electronic address: sroy@iopb.res.in magnetic susceptibility and correlation length of the system, t 

^Electronic address: somen@iopb.res.in is the temperature-like variable with the critical point at t = 0. 
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the correlation volume, yielding, 

[(6T C )% V ~ N~ 2 ' p (3) 

An immediate consequence of this is that R x of Eq. 
approaches a constant as N — > oo indicating complete 
absence of self-averaging at the critical point in a ran- 
dom system. These predictions have been verified for 
various types of relevant and irrelevant disorders and 
also with canonical (ensemble of fixed concentration of 
disorder) and grand canonical disorder at the random 
critical point [111110 . 



II. SMALL- WORLD NETWORK 

Over the last few years, small-world networks (SWN) 
have emerged as a new class of graphs with character- 
istic statistical properties defined over the ensemble of 
networks. Starting from an Euclidean lattice, one may 
obtain an SWN by rewiring the original lattice or by 
the addition of random long range bonds, even with a 
sufficiently small probability pjfj H @- The network 
is so named because any two points far away on the 
underlying lattice can be bridged, on the average, by 
a finite number of connections. It has been observed 
that for an Ising model defined on such an ensemble of 
graphs, this set of random bonds changes the critical be- 
haviour to the mean field type in all the models studied 
so far dEl El. This, then, by definition, makes 
this set of random bonds, added to the underlying "pure" 
lattice, a "relevant" variable. There is albeit a debate on 
the crossover exponent for p at the pure critical point, in 
the limit p^0[ll|,ll2j. Since quenched averaging is im- 
portant, replica trick has been resorted to [10. Over- 
laps and other quantities of interest for gene ral quenched 
random systems have also been studied |14j. 

The aforementioned result of non-self-averaging would 
imply a strong influence of the network i.e. a network to 
network variation of an extensive quantity at the shifted 
critical point. This aspect of disorder of the small- world 
networks is the primary motivation of the study reported 
here. 

With this background we set to check the behaviour 
of Rx for various X for SWN. Some of the major dif- 
ferences with respect to the previous studies of random 
systems may be mentioned here. By construction, the 
random bonds in SWN introduce long-range interactions 
unlike the short range models studied so far. The self- 
averaging behaviour of long range cases is of importance 
0,0/ because of the specialities known, e.g., for disorder 
with long range correlations lla l. We would also like to 
note that compared to many other random systems, the 
exact values of the changed exponents are known in the 
SWN case. It is well established that the shifted critical 
behaviour has a — and therefore D = 2. 



III. THE ISING MODEL ON A SMALL- WORLD 
NETWORK 

Let us think of the ferromagnetic Ising model with 
spins Si = ±1 at each site i of a network based on an 
Euclidean lattice of N points similar to Ref. [1 IHl 
Il2l IT^ . IbH ITt| . Its Hamiltonian is taken as 

h({s}) = -jJ2^j- jJ2 s ^> ( 4 ) 

where (ij) are the nearest-neighbours on the lattice , (ij) 
are the long distance neighbours along the random bonds 
of {S} added for the network and J > 0. The Hamil- 
tonian, being dependent on the set {S}, is random for 
a given configuration of the spins {si}. Any physical 
property X of the model, therefore, requires an overall 
averaging over the ensemble of networks. For our stud- 
ies we need the shifted critical point and the sample to 
sample fluctuations of various quantities at this critical 
point. Since the random bonds are known to be relevant, 
for all d < 4, there is no loss of generality if the SWN is 
built from a one dimensional lattice. 

We start with the Ising model on a SWN in ID. In 
our model each site on the lattice with an Ising spin has 
random links to two distant spins such that no two spins 
are connected by more than one link. All links of equal 
strength. Thus we have a "canonical" scenario since the 
number of links at each site is fixed. Hence no extra 
normalisation factor is needed in the long range part of 
the Hamiltonian of Eq0J In the present work we chose 
J/ks = 1 where ks is the Boltzmann constant. 



IV. SIMULATION DETAILS 
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[T/Tc - 1 ] N 1 /2 

FIG. 1: Plot of the data-collapse of the Binder Cumulant 
versus scaled temperature for various TV 

Data were taken at T = 2.85 (close to the estimated 
critical temperature) and Xi C and the Binder cumulant 
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U m = [< m 4 >/< m' 



2i 
> av 



3 were calculated, using 



the single histogram reweighting technique |l7j . We ex- 
amined lattice sizes N = 100, 500, 1000, 2000, 3000 in our 
Monte Carlo simulations. We studied 1535 samples for 
N = 100 to 517 samples for N = 3000 using 10 3 equili- 
bration and 10 6 MC steps for each TV. Data were taken 
at intervals of 10 3 MC steps . 

A data-collapse of U m with finite size scaling variable 
N X /"(T-T^ /T c would give T c (the infinite lattice critical 
temperature) and v. By using the data-collapse method 
of Ref. HI we obtained T c = 2.83(2) and 1/u = 0.50(1) 
The value of v is consistent with previous results [lJl • 
The resulting collapse is shown in Fig. We also investi- 
gated similar plots for \ and C after averaging over many 
realisations of disorder (not shown here), with same re- 
sults. 

Further proof of the mean-field nature of the transition 
comes from the comparison of the data with the mean- 
field form of U m ■ To evaluate the mean field form of U m 
we use the mean-field form of the magnetisation per spin 
m probability distribution in the critical region |19| 



Pjv(to) oc exp [— Niait^m 2 + a2m 4 )] 



(5) 



with t 7v being the critical temperature of a lattice of 
size N and oi, a 2 being constants. By replacing fa — 
(a2^V) 1 / 4 m, we find U m where the averages are obtained 
by integrating rh from — oo to +oo with the weight 



P/v(m) oc exp [— bi(x — b^rh 2 — m 4 



(6) 



where b\ = ai/a2 1 ^ 2 , x is the finite size scaling variable 
[(T-iy/TjiVVP with bib 2 taking care of the finite size 
shift of the critical temperature. The solid curve in Fig. 
^is obtained with b\ = 1.7, bi = 0.72. 

We find good data collapse by using t = T — T c even 
though finite size scaling is supposedly better with the 
use of t = T - TAN) after finding out the T C (N) for 
every sample size[2(|. This is because our system, as we 
show, is self-averaging and this method should be more 
pertinent for non- or weakly self-averaging systems. To 
investigate the distribution of pseudo-critical tempera- 
ture, (3 c {i, N) (the temperature at which the specific heat 
of sample i of size N is a maximum), data were taken 
at T = 2.85 and (3 c {i,N) for various N were calculated 
using the histogram method 0. The distribution of 
(3 c (i,N) for N = 10,50,100 is constructed.We studied 
3291 lattice samples for N = 10, 1645 lattice samples for 
N = 50 and 1535 lattice samples for N = 100. We find 
that the inverse critical temperature (3 c (i,N), scales as 



Kf3 c (i,N)-f3 c (N)) 2 



N 



-2/v 



(7) 



which is consistent with Eq. but one also needs 

a scale factor TV" 1 / 2 for the probability distribution 
P(f3 c (i, N)). Fig. [5] shows the data collapse of this 
distribution. The data-collapse is best achieved with 
v = 2.00(3) which is consistent with the value of v ob- 
tained in the data collapse shown in fig. ^ We then 
studied R m , R x and Vc at r c °° for the above lattice 
sizes. About 56440 lattice samples for N=10 to 1000 
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FIG. 2: Plot of the distribution of pseudocritical inverse tem- 
peratures flc(i,N) for various TV 





FIG. 3: (a) R m , R x versus N at T c and (b) Vc versus N at 
Tc-The straight lines show straight line fits to R m and R x . 



lattice samples for N=3000 were studied. For each sam- 
ple we used 10 3 equilibration and 10 5 MC steps. Data 
were taken at intervals of 10 3 MC steps. The data is 
fitted to the form Rx = AxN px where Rx is the rel- 
ative variance for m and \- The values obtained are 
Pm = —0.96(9), and p x = —0.94(8). Thus x an d m are 
strongly self-averaging . The singular part of energy can- 
not be filtered out and hence the behaviour of Ve can 
not be predicted decisively. We see in Fig. 3 that Vc 
is a constant as expected and hence C is also strongly 
self- averaging . 



V. DISCUSSION 

The fact that the peak scales inversely as the width 
shows that despite the fluctuation in pseudo-critical tem- 
peratures, the distribution approaches a S— function. As 
a result the critical temperature of a large iV network 
can be thought of as the average of pseudo-critical tern- 
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FIG. 4: Data collapse of C and \ on considering a single 
realisation of randomness 



peratures of the small sub-networks. This averaging out 
is tantamount to self-averaging .This is in marked con- 
trast to the other cases of random systems studied so 
far 0, 0, 0, H • It is tempting to conclude that in ad- 
dition to relevant randomness, a broad distribution of 
pseudo-critical temperature is a requirement for non-self- 
averaging. 

Whilst in the present work we have used a "canoni- 
cal" ensemble with a fixed number of bonds, small scale 
simulations of the Ising model on a SWN in a "grand 
canonical" ensemble, where the number of bonds can 
vary, also indicated self-averaging as found here. 

In case of a strongly self-averaging system, a typical 
sample should be a representative of the average. We 
observe good data collapse with even a single realisation 



of disorder (as shown in Fig. Thus, in the light of the 
present work, in such situations, an annealed averaging 
as done in ref 0] should work well. Consequently no 
extra order parameter like approach should be needed 
for networks . 

It is not clear if this feature of strong self-averaging is 
a consequence of a — 0, in which case it could be true for 
all relevant disorder problems with mean-field behaviour. 
An extension of the RG argument [l| to encompass situ- 
ations with sharp limit of T c distribution and long range 
interactions, may shed light on this. Whether this re- 
sult on the disorder aspect of a network is important 
in other real life situations like the railway network |21j 
needs further study. 

To conclude, we investigated the self-averaging be- 
haviour of the Ising model on a small world net- 
work. The distribution of (3 c (i,N) is found to become 
sharper as N — > oo with the fluctuation decaying as 
[(<%(i, A)) 2 ] av ~ N- 2 ^. The data collapse of various 
physical quantities both for a single realisation of disor- 
der and after averaging over many disorder realisations 
showed no significant difference. At T c °°, the relative 
fluctuations R rn , R x for magnetization and susceptibil- 
ity are found to behave as R m ,R x ~ A -1 while the 
variance Vc for the specific heat approaches a constant 
for large N. Hence the system is strongly self-averaging 
in the critical region in spite of relevant randomness. 
Our results have the following implications. From the 
small-world networks perspective, the random or statis- 
tical features of the network do not play a role in the 
long range behaviour at the critical point so that one 
may replace the ensemble of small- world networks by a 
single average network. 
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